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We analyze the dynamics of a nano-mechanical resonator coupled to a single-electron transistor 
(SET) in the regime where the resonator behaves classically. A master equation is derived describing 
the dynamics of the coupled system which is then used to obtain equations of motion for the 
average charge state of the SET and the average position of the resonator. We show that the 
action of the SET on the resonator is very similar to that of a thermal bath, as it leads to a 
steady-state probability-distribution for the resonator which can be described by mean values of 
the resonator position, a renormalized frequency, an effective temperature and an intrinsic damping 
constant. Including the effects of extrinsic damping and finite temperature, we find that there remain 
experimentally accessible regimes where the intrinsic damping of the resonator still dominates its 
behavior. We also obtain the average current through the SET as a function of the coupling to the 
CZ3 ' resonator. 



I. INTRODUCTION 



Recent developments in fabrication have made it possible to produce nano-electromechanical systems in which 
mechanical resonators with frequencies in the range 100 MHz-1 GHz are coupled electrostatically to mesoscopic 
conductors ^^^^ The dynamics of nano-electromechanical systems have attracted considerable interest because 
of the novel transport mechanisms they can give rise tof&SilfiiAi their interesting non-linear and chaotic properties^iSii^ 
and their applications as fast and ultra-sensitive sensorsi^ii Furthermore, when cooled to ultra-low temperatures, high- 
frequency resonators are expected to display quantum mechanical behavior and so nano-electromechanical devices 
I ' could be used to explore the cross-over from quantum to classical behavior in mechanical systems^i^i^iii 

TheoreticaljiSiiSiSSi and experimental^ work has suggested that a single-electron transistor (SET) could form the 
basis of an extremely sensitive motion detector for nano-mechanical resonators. Therefore, the question of what the 
£NJ ' potential sensitivity of such a displacement detector will be, and how it is determined by the underlying dynamics of 
W) • the coupled SET-resonator system, is of practical as well as theoretical interest. 

When a mechanical resonator, coated with a thin metallic layer, is placed next to a SET, it acts like a capacitor 
whose presence affects the current flowing through the SET. The capacitance, and therefore the current, depend on 
the position and hence the motion of the resonator. The sensitivity with which the displacement of a resonator can be 
measured from the current through the SET is limited by two sources of noise. Firstly, the shot noise in the current 
through the SET limits its sensitivity. Secondly, the fluctuating charge on the SET island will act as a fluctuating 
force on the resonator, resulting ultimately in displacement noise. Thus far, the sensitivity of the displacement 
I ■ detection which can be achieved has been gauged by treating the resonator and the SET as two distinct elements^ 
Although this approach provides a reasonably accurate estimate of the sensitivity, it neglects correlations between the 
displacement noise of the resonator and the shot noise in the current. In fact, the coupled SET and resonator form 
an intrinsically interesting dynamical system in their own right. By studying the full coupled dynamics of the system 
a broad understanding of the underlying physics can be achieved which can then be applied to obtain important 
insights into when the system would function most sensitively as a measuring device. 

In this article, we analyze the coupled dynamics of the SET-resonator system, treating the resonator as a classical 
harmonic oscillator. In particular, we examine the effect on the motion of the resonator of electrons tunnelling through 
the SET and the concomitant effects of the resonator motion on the average current of the SET. More generally, our 
analysis of the coupled dynamics provides useful insights into the dynamical behavior of nano-electromechanical 
systems. 

The dynamics of a quantum harmonic oscillator coupled to an electrical tunnel junction was studied recently by 
Mozyrsky and Martini^ and by Smirnov et al. 22 In the regime where the drain-source voltage across the junction, 
Vds, is much larger than the energy quanta of the oscillator, Hloq, Mozyrsky and Martin found that the dynamics 
of the density matrix of the oscillator reduced to a form very similar to that of the Caldeira-Leggett models for 
an oscillator in contact with a thermal bath. From the connection with the Caldeira-Leggett model, Mozyrsky and 
Martin were able to show that the tunnelling electrons caused rapid dephasing of the oscillator, acting like a thermal 
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bath with an effective temperature and associated damping constant. 

For a classical resonator coupled only to a SET, we find that the electrons affect the resonator in a similar way to 
the quantum system considered by Mozyrsky and Martini For sufficiently weak coupling between the resonator and 
the SET, the electrons act like a thermal bath which heats the resonator to an effective temperature. However, if the 
coupling between the SET and the resonator is increased the behavior of the resonator begins to deviate significantly 
from that of an oscillator in contact with a thermal bath. Including the effects of extrinsic damping of the resonator, 
we find that the effective temperature of the resonator is suppressed. 

The classical regime for the resonator which we consider here will be the appropriate description for systems where 
eVds S> Hujo and the model we use to treat the electron tunnelling processes in the SET will require eVds ^> ksT e ^ 
where T e is the electron temperature (assumed to be the same in the SET island and leads). As the electrons passing 
through the SET heat the resonator to an effective temperature which is proportional to eV^s m the absence of extrinsic 
damping, this means that the resonator effective temperature can be much larger than the electron temperature. 
We will assume possible overheating of the SET can be neglected so that this temperature imbalance can be 
maintained and the above model assumptions continue to hold. 

The outline of this paper is as follows. In Sec. II, we review the properties of the SET with a fixed voltage gate 
and the master equation method which can be used to determine the average current. Then in Sec. Ill, we generalize 
the results for a fixed voltage gate to the case where the gate is a nano-mechanical resonator. Initially, we analyze 
the dynamics of the resonator coupled to the SET ignoring external influences. We derive the equations of motion 
for the probability distribution of the system and show that, with certain approximations, the average properties 
of the system can be described by simple, physically intuitive equations of motion. We investigate the steady-state 
properties of the probability distribution for the coupled system by numerical integration of the master equations that 
govern the dynamics. 

In Sec. IV, we consider the effect of damping in the resonator, arising from sources other than the interaction with 
the SET, and the effect of finite background temperature. We examine the current characteristics of the system in 
Sec. V, and obtain the average current through the SET when it is coupled to the resonator. Then in Sec. VI, we 
discuss how our results could be investigated experimentally. Finally, in Sec. VII, we draw our conclusions and discuss 
the more general implications of our results for understanding the dynamics of nano-electromechanical systems. 

In the appendices we give further details about aspects of our calculations. In appendix A, we derive the Hamil- 
tonians and electron tunnelling rates for the coupled SET-resonator system. We give details of the method used to 
integrate the master equations for the SET-resonator probability distribution in appendix B. The effect of higher 
vibrational modes of the nano-mechanical resonator is discussed in appendix C. Finally, in appendix D we estimate 
the typical experimentally achievable values of the parameters appearing in our model. 

II. CHARGE DYNAMICS OF A SET COUPLED TO A FIXED GATE 

The development of the single-electron transistor (SET)2£ and its subsequent refinement offers the prospect of 
carrying out electrometry with unprecedented accuracy^ A SET consists of a small metal island joined to leads 
by two tunnel junctions. The small capacitance of the metal island can give rise to a Coulomb blockade where the 
charging energy of a single electron dominates over thermal effects at low temperatures and the only accessible state 
of the system has a fixed number of electrons on the island, N. An external voltage-gate can be used to induce a 
polarization charge on the metal island which makes the states with N and N + 1 electrons degenerate, and thus 
allows current to flow through the device. In the low temperature regime, states with more than N + 1 or fewer than 
N electrons have much higher energies because of charging effects and so play no role in the dynamics of the system.— 

Before we go on to consider the coupled dynamics of the SET-resonator system, we briefly review the properties 
of the SET when coupled to a fixed voltage-gate. The circuit diagram for a nano-mechanical resonator coupled to 
the island of a single-electron transistor is shown schematically in Fig. 1. Neglecting the effect of the dynamics, the 
resonator acts as a simple voltage-gate with a fixed position. Indeed, the dynamics of the resonator, which we consider 
in the next section, generally acts as a small perturbation on top of the effect of a static voltage gate. 

We simplify the analysis a little by assuming that the SET is symmetrical, as shown in Fig. 1. The two junctions 
linking the metal island to the leads have the same capacitances, Cj, resistances, R, and a drain-source voltage Vds 
(which we take to be positive throughout) is applied across the island symmetrically The static resonator acts as the 
gate of a capacitor, C g , with a voltage V g applied to it. The gate voltage can have a strong influence on the charge 
dynamics of the electrons on the island because of the polarization charge it induces. 

The energy of N charges on the SET island with a fixed gate is given, up to an unimportant constant factor, by^i 



H N = ^-{N 2 -2NN g ), 



(1) 
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FIG. 1: Circuit diagram for the coupled SET-resonator system. The SET island lies between the two tunnel junctions which 
are assumed to have the same capacitance Cj. Under the drain-source voltage shown, the electrons tend to flow through the 
SET island from right to left. 



where Cs = 2Cj + C g is the total capacitance of the SET, and N g = —C g V g /e is the polarization charge induced by 
the gate. 

In the regime where e 2 / (2Cs) S> A:bT £ , and where the bias voltage is not too large, the electron tunnelling processes 
in the SET are restricted to those between N and N + 1 electrons on the islandii2i22i2I The tunnelling dynamics of 
the electrons can be modelled by a simple master equation known as the orthodox models This approach is valid 
provided that the effective tunnel junction resistances satisfy R > h/e 2 , the charging energy dominates over the 
thermal energy, and the SET is biased well away from the Coulomb blockade regions so that co-tunnelling can be 
neglected. 

We now outline the form of the master equation and the way in which the steady-state current is obtained for the 
SET with a fixed gate. The analysis of the coupled dynamics of the SET-resonator system is then carried out by 
generalizing these results to include the motion of the resonator. 



A. Master Equation 

Within the orthodox model^i the state of an ensemble of SETs for which only the two charge states N and N + 1 
are accessible can be described by the probabilities PN(t) and PN+i(t), of finding N or N + 1 electrons on the island 
at time t, respectively. These probabilities evolve in time according to a pair of coupled master equations, 



dP, 



N 



dt 

dPN+i 



= -(ti + r+)p N + (v+ + r R )p N+1 (2) 
iU = (r- L +T+)p N -(r+ + T R )p N+u (3) 

where are the tunnelling rates forwards (+) or backwards (— ) across the two junctions L and R between the 

island and the leads, as illustrated in Fig. Here we use the convention that the forward direction is from right to 
left, matching the direction set by the drain-source voltage. The tunnelling rates are calculated using the orthodox 
model and at finite temperature, T e , are given by^ 

where AE^ R s is the initial free energy (of the SET and leads) before tunnelling through the L(R) junction, minus the 
final free energy after tunnelling forwards (+) or backwards (— ). At T e — 0, this expression reduces to the simplified 
form, 

r± = ^R- AE L(R) Q ( AE i(R))> ( 5 ) 
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where Q(x) is the Heaviside step-function. 



B. SET Current 



The current through the SET depends on the tunnelling rates through the two junctions. In the steady-state, 
current conservation means that we need only ever consider electron tunnelling processes across one junction. When 
the Coulomb- blockade limits the accessible charge states of the SET island to N and N + 1, the steady-state current 
can be written as^ 

(7) = e {T+P N+1 - T- L P N ) (6) 
= e (T+P N - T-P N+1 ) , (7) 

where (• • •) denotes an ensemble average and where P/v and Pn+i = 1 — -P/v are the probabilities of having N and 
N + 1 electrons on the island, respectively. Substituting in the correct expressions for the tunnelling rates, we find 

(I) = -P N P N+ll (8) 

n 

where r t = Re/V^g is a measure of the average time between tunnelling events. 



III. DYNAMICS OF THE COUPLED SET AND RESONATOR 



Unlike a fixed voltage-gate, a movable voltage-gate, such as a nano-mechanical resonator, responds dynamically to 
the force exerted on it by electrons tunnelling through the SET island. Furthermore, the motion of the resonator acts 
back on the SET by altering the tunnelling rates of the electrons. The coupled dynamics of the SET electrons and 
the resonator eventually results in a steady-state, even without including the effects of external sources of dissipation 
on the resonator. 

The dynamics of the SET-resonator system are described by generalizing the master equation for the SET charges 
to include the position and velocity of the resonator. Crucially, the tunnelling rates for electrons in the master 
equation must be re-calculated to include the effect of the resonator's motion. Although the master equation for the 
SET-resonator system is apparently much more complicated than that for the fixed-gate case, the dynamics of the 
average resonator position and SET charge are readily determined once certain simplifying approximations are made. 
Furthermore, the master equation can be solved numerically and it is found that for sufficiently weak coupling the 
steady-state probability distributions are approximately Gaussian in form. The effect of the SET on the motion of 
the resonator can then be gauged very accurately from the first and second moments of the steady-state probability 
distribution. 

We treat the nano-mechanical resonator as a single-mode harmonic oscillator, as it turns out that the interaction 
between the SET and the fundamental flexural mode of the resonator dominates the behavior. However, if we assume 
that interactions between the various vibrational modes of the nano-mechanical resonator can be neglected, then the 
results we obtain for a single mode can be generalized to include higher frequency modes. 



A. Master Equation 



The probability distribution for the SET-resonator system is a function of the oscillator position, x, and velocity, 
u, as well as the charge state of the SET island. We can characterize the state of the system by the probability 
distribution pair, P^{x,u]t) and PN+i(x,u;t), which give the probabilities of the oscillator being at position x, with 
velocity u, and having N or N + 1 electrons on the island, respectively, at time t. At T e = 0, the equations of motion 
for these probability distributions take the form, 

dP N (x,u;t) 

= {H N ,P N {x,u;t)} (9) 



dt 



~ [e(AE+)AE+ + Q(AE R )AE n ] P N+1 (x, u; t) 
- 1 ^[e(AE+)AE+ + Q(AEZ)AEl]P N (x,u;t) 
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dP N+1 (x,u;t) 

= {H N+u P N+1 (x,u;t)} (10) 

at 

+ ^[e(AE+)AE+ + Q(AE-)AEz]P N (x,u;t) 

[e(AE+)AE+ + Q(AE R )AE R ] P N+1 (x, u; t), 

where Hn(n+i) is the Hamiltonian for the SET-resonator system with N(N + 1) electrons on the island, and {•, •} 
is a Poisson bracket. The quantities AE^ R ^ are the total free energy differences for an electron tunnelling forwards 
(+) or backwards (— ) across the left-hand (right-hand) junction. 

The Hamiltonians for the coupled SET-resonator system with N and N + 1 electrons on the SET island are readily 
obtained treating the resonator as a single- mode harmonic oscillator, with angular frequency ujq, centered a distance 
d away from the SET island. Details of the calculation are given in appendix A. The Hamiltonians take the form 

H N = E c SN+^- + ^mtu 2 x 2 (11) 

p 2 1 1 
Hn+i = -E C SN + [- -muj 2 (x - x ) 2 ~ -muj 2 xl (12) 

£uFYh — £ 

-muj 2 xl(N g -SN- 1/2), (13) 

where SN = N g — N — 1/2 and xq = —2E c N g /(mujQd). The length-scale Xq has a simple physical interpretation: it 
is the distance between the equilibrium positions for the resonator with N and N + 1 electrons on the SET island. 
The energy differences which determine the tunnelling rates are given by 

AE+ = -AE- = -2E c 5N -mufcox-mujfaliNg- 5N -1/2) + ?-^- (14) 

AE+ = -AE R =2E c 6N + muj 2 ,xox + muj 2 xl(Ng- <«V - 1/2) + (15) 

as derived in appendix A. Notice that at T e = 0, for a given displacement of the resonator and voltage bias, tunnelling 
across each junction will only be able to occur in one direction. The equations of motion for the probability distribu- 
tions can be extended to finite temperatures by replacing the tunnelling rates by the appropriate finite temperature 
form, as given in Eq. ijlf. 



B. Mean-Coordinate Dynamics 

Integration of the coupled master equations gives a complete description of the ensemble dynamics. However, the 
equations of motion for the probability distribution are complicated and it is not easy to extract an overall picture of 
the dynamics from them. Nevertheless, by making a few physically motivated approximations, it is possible to obtain 
equations of motion for all of the first and second moments of the distributions. 

In order to simplify the analysis, we assume that SN, Vd s , and the range of x for which the probabilities, Pn{x, u] t) 
and Pn+i(x, u; t), are non-negligible are such that AE^,AE R are always positive. Such an approximation is valid 
so long as the drain-source voltage is the dominant energy scale so that eVds ^> itlujqX 2 ,, EcSN and the distributions 
Pn(x, u; t) and Pn+i(x, it; t) are strongly peaked at x ~ x$. 

Taking into account this approximation, we need only consider two out of the four tunnelling processes, as the step 
function constraints will always be satisfied for the processes in the '+' direction, and never satisfied for processes in 
the '— ' direction. Thus, we obtain a simplified pair of coupled master equations, 

dP N 2 dP N dP N 1 , 2 p \ mrs 

' —2 (E L P N+ i - E R P N - muj xoxP) (16) 



u- 



dt du dx Re 

dPN+l 2( n<9P/V+1 dP N+1 1 , 2 u \ fT7\ 

= u>$(x- xo)—q^ u— — — (E L P N+1 - E R P N - rmj£x xP) , (17) 

where P(x, u; t) — Pn(x, u; t) + P N+ i(x, u; t), and we have defined E L = —2E C 5N — muj 2 xl(Ng — SN — 1/2) + eVds/2 
and En = 2E C SN + muj 2 XQ(N g — SN — 1/2) + eVd s /2 (i.e., El and Er are the position-independent parts of AE£ 
and AE R , respectively). 

Adding together Eqs. i|ltj|) and l|17|) . the tunnelling terms cancel out and we obtain the simplified equation of motion 
for the full probability distribution. 

dP 2 dP OP dP N+1 
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We can now obtain equations which describe the dynamics of the ensemble-averaged position and velocity by multi- 
plying Eq. (|18|) by x and u respectively and then integrating over both x and u, 

(u(t)) = -u>l(x(t)) +ujlx (P) N +i(t) (19) 
(*(*)} = («(*)> (20) 

where the averages are defined by 

(•••) = J dx J du{-- -)(P N+1 (x,u;t) + P N (x,u;t)) 

and 

(P)n+i = J dx J du P N+1 (x,u;t). 

We can obtain a closed set of equations by integrating Eq. (jl7[) over x and u to obtain, 

Re 2 {P) N+1 (t) =E R - eV ds (P) N+1 (t) - mu 2 x {x(t)). (21) 



1. Fixed Point Analysis 

The steady-state behavior of the system is described by the fixed-points of the mean-coordinate equations 119I20I21|1 . 
Setting (u) — (x) — (P)n+i = 0' we nn< ^ that there is a single solution at 

(u) ip = (22) 
{x)fp = x {(P) N+1 ) {p (23) 
2E C SN + mu$x%(N g - 8N - 1/2) + eV ds /2 
eVds - muj^xl 



((^iV+ljfp = — — , 2 2 • I 24 ) 



In order to make the physical content of the mean-coordinate equations H19I20I21(I more apparent, we transform 
them to variables centered on their fixed points, 

(x) = (x)-(x) fp (25) 
(P) N+1 = (P) N+1 -((P) N+1 ) ip , (26) 

and we also re-express them in terms of the dimensionless frequency and coupling strength parameters e = LUQT t and 
k = mujQxl / (eVds) , respectively. Hence we find, 

(P) N+1 - —{£) - (P)jv+i (27) 
x 

($) = (u) (28) 

(u) = e 2 (x (P) N+1 - (x)). (29) 

where the time coordinate has been scaled by the tunnelling time, r t . 

The scaled parameter e compares the relative time-scales of the oscillator period and electron tunnelling time. 
Generally, the oscillator motion will be much slower than the electron tunnelling time and so we expect e <C 1 for 
most practical situations. The parameter k is a measure of the interaction strength between the SET and the resonator 
in terms of the energy-scale defined by the drain-source voltage. In practice, the coupling between the resonator and 
the SET which can be achieved will be weak so that n <C 1, as discussed in Sec. VI. 

The mean-coordinate dynamics about the fixed point is given by a set of three coupled, linear, first-order differential 
equations, and the stability of the fixed point can be determined from the behavior of the associated eigenvalues. 
Although the eigenvalues of a system of three equations can be determined analytically, the expressions obtained 
are algebraically rather complicated. However, the stability of the system is determined by whether or not the real 
part of any of the eigenvalues becomes positive, in which case the system can become unstable. A simple numerical 
calculation shows that the SET-resonator system is entirely stable for k < 1, but for k > 1 it seems that our equations 
predict that the system becomes unstable. However, any potential instability in the motion of the resonator would 
clearly involve processes which lie outside the domain of our approximate description and so we cannot draw any 
conclusions about the overall stability of the system within the existing model. 
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FIG. 2: Comparison of the approximate analytical equations of motion for the resonator with the exact equations, integrated 
numerically. In this plot k = 0.1 and e = 0.3, so the approximate equations are close to the limit on their validity (e <C 1). 



2. Analytical Description 



The set of linear equations describing the mean-coordinate dynamics of the SET-resonator system can be solved 
analytically, and the general solution obtained. However, as with the eigenvalues, the general solution is algebraically 
too complicated to be useful in developing an intuitive picture of the dynamics. However, in the regime where e<l 
and K<la very good approximation to the full solution can be derived as an expansion in e. 

The approximate solutions for the trajectory of the system take the form, 



(P)N+l(t) 



-Ke 2 t/2r t 



x(0) cos [yl — K et/rt) + 



u(0) 



y/1 — K € 



sin (vl — k et/rt) 



e - K e 2 i/2r t _^ 
X 



x(0) cos - K et/rt) + U ^ sin (Vl - K et/r t ) 

\/l — k e 



(30) 



(31) 



with the initial conditions (x(t = 0)) = i(0), (u(t = 0)) = u(0) and (Pjv+i)(i = 0) = Pn+i{0). The approximate 
trajectory is in fact very close to the full solution, as is shown in Fig. [5] for k = 0.1 and e = 0.3. The value of e in 
the plot is at the limit of the range of validity of the approximate equations, but the agreement remains good. For 
values of e < 0.1 the approximate and exact curves overlay each other so well that the two curves cannot easily be 
distinguished. 

The approximate analytical solution for the dynamics of the oscillator given by Eq. (|30|) can be obtained from 
the equation of motion for a damped harmonic oscillator, with damping constant ne 2 /rt, and (angular) frequency 
VI — K e/r t . The shift in frequency is due to the interaction between the oscillator and the tunnelling electrons. 23 
We can think of 7, = ne 2 /rt as an intrinsic damping rate for the SET-resonator system, to be distinguished from 
extrinsic damping arising from sources other than interaction with the SET (see subsection II V A|l . 



C. Numerical Solution of the Coupled Master Equations 

As well as extracting mean-coordinate properties, it is possible to integrate the equations for the probability 
distributions [Eqs. (|16f) and l|17l) ] numerically, using a method which is outlined in appendix B. Numerical integration 
allows us to demonstrate that the coupled master-equations do indeed lead to a steady-state probability distribution 
and also to explore how this steady-state is reached. 

Figure |2| shows the dynamics of the variance, Sx 2 , of the full probability distribution, P{x,u\t). The initial 
probability distributions were chosen to be sharply peaked in phase space. For simplicity we have performed all 
the calculations at a gate-voltage such that El = En = eVd s /2, so that the N and N + 1 charge states would be 
degenerate in the limit K = 0. After some initial oscillations, Sx 2 settles down to a steady value. The rate at which 
the steady-state is achieved is ji, just as we would expect for a damped harmonic oscillator. Notice, however, that all 
the lengths have been scaled by xq, which itself depends on k. Thus, the absolute magnitude of the variances increase 
with k in contrast to the scaled values which actually decrease. 

The steady-state probability distributions P(x), for e — 0.3 and k — 0.1,0.05, and 0.01, respectively, are shown in 
comparison with the appropriate Gaussian fits in Fig. 0] The steady-state distributions are indeed strongly peaked, 
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600 800 



600 800 



i/x, th, 

FIG. 3: Behavior of the variances, 5a: 2 and 8u 2 , as the steady-state is approached. In panels (a) and (c) k = 0.05, and in panels 
(b) and (d) k = 0.1, with e = 0.3 throughout. The dashed lines represent fits to an exponential time-dependence obtained using 
the initial variances, the steady-state values calculated using Eqs. 1121 and l|43|l . and the appropriate decay rate, ji — e 2 /t, in 
each case. 

as we assumed earlier, and the overall distribution is very nearly Gaussian as can be seen from the plots which are 
generated from the average position and associated variance in each case. The deviation from a Gaussian distribution 
is most apparent for k = 0.1 and the agreement improves as the coupling is reduced. The sub-distributions Pn{x) and 
Pn+i(x) for k = 0.1 are shown in Fig. ^d). They are clearly centered on different points and also closely resemble 
Gaussians individually. However, since one cannot add two displaced Gaussians to obtain a single Gaussian curve, it 
is clear that the curves are not exactly Gaussian in form. 

D. Variances of the distribution 

Since the steady-state distribution is close to Gaussian for k <C 1, the state of the resonator in this regime will 
be almost completely specified by the average values of its coordinates and their variances. Making use of the fact 
that there is indeed a well-defined steady-state, we can derive analytical expressions for the variances in the resonator 
position and velocity. 

Multiplying Eqs. 116ll7f) by the appropriate quantity and integrating over the position and velocity of the resonator, 



we obtain the following closed set of equations of motion: 

(xu) = -ul{x 2 ) + (u 2 ) +w%x (x) N+1 (32) 

(x 2 ) = 2{xu) (33) 

(u 2 ) = —2u) (xu) + 2x u1q(u) n + x (34) 

(x) N = (u) N + [El(x)n+i - E R (x) N - mw 2 i (i 2 )] (35) 

(x) N+1 = (u)n+i - 77-7 [E l (x) n +i - E R (x) N - mw 2 i (a; 2 )] (36) 

(u) N = ~lu 2 (x) n + [E l (u) n +i - E R (u) N - mu! 2 x (ux}] (37) 

(u)n+i = -Vo(x)n+i + ll>Ix (P) n+1 - — [E L (u) N+1 - E R (u) N - muj 2 xo(ux)] , (38) 



where (• • -)jv and (• • -)n+i imply averages over the sub-distributions Pn(x,u) and Pn+i(x,u), respectively. 

In the steady-state, all time derivatives will be zero and hence we can infer the following properties of the stationary 
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FIG. 4: Steady-state probability distributions as a function of position. The full probability distribution, P(x), centered on the 
fixed point xo{P) jv+i, is shown in (a) with k = 0.1, (b) with k = 0.05, and (c) with n — 0.01, where e = 0.3 in each case. Each 
of the numerical curves is compared with a Gaussian fitted using the fixed-point position and the position variance, obtained 
from Eqs. (I23H and 14211 respectively. In (d) the sub-distributions Pn{x) and Pm+i{x) for k = 0.1 and e = 0.3 are shown. 



distribution 



{x)n+i = x Q (P) N+ i (39) 
(u 2 ) - loq(x 2 ) = ~uj 2 x (x) n+ i (40) 

(x 2 ) = ^(P)n+i, (41) 

and (xu) — (u)n — (x)n = (u)n+i = 0. Taken together with the steady-state solutions of the mean-coordinate 
equations, (x) — xq(P)n+i and (u) — 0, we obtain expressions for the variances in position and velocity for the 
stationary probability distribution. Thus we find, 

fe 2 = ?0 {P)n{P)n+i = ^1 {P)n{P)n+1 (42) 

Su 2 = ^(i- K ){P) N {P) N+1 =*V±(i- K ){P) N {P) N+1 . (43) 
k m 

The accuracy of these expressions, based on the assumption that only forward electron tunnelling occurs, can be 
tested by comparing them with the results obtained from numerical integration of the master equation. Fig. 
compares the analytical result, Eq. 142|) . with the numerical result as a function of k. It is clear that there is excellent 
agreement for k < 0.1, but for larger values of k the validity of the analytical expression breaks down. 

In order to understand the expressions for the variances, Eqs. Q42J) and (|43(l . we can compare them with those for 
a damped harmonic resonator in contact with a heat bathiS For an oscillator in contact with a bath at temperature 
T, we have equipartition of energy: muS^Sx 2 — m5u 2 — k^T. Comparison with Eqs. I|42|l and l|43|) shows that for 
k <C 1 the steady-state probability distribution for the resonator is very close to a thermal one, and we can identify 
an effective temperature, T e g, by analogy with the thermal case, 

k B T cff = eV ds (P) N (P) N+1 . (44) 

The electrons passing through the SET give rise to a fluctuating force on the oscillator as its equilibrium position is 
shifted back and forth by xq, but because the oscillator acts back on the motion of the electrons through the tunnelling 
rate, the motion of the oscillator is also damped and a steady-state is achieved. For k« 1. the steady-state is very 
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FIG. 5: Comparison of the variance, 6x 2 , calculated using Eq. 1421 (curve), and the results of a numerical integration (points). 
Both calculations are performed for El = Er. 

close to a thermal state with an effective temperature which depends on properties of the SET alone: the drain-source 
voltage applied across it and the average probability of finding N or N + 1 electrons on the island. In particular, the 
effective temperature docs not seem to depend at all on the strength of the coupling between the resonator and the 
SET. In fact, the strength of the coupling k does affect the effective temperature, as the values of (P)n+i and (P)n 
depend (albeit weakly) on k for n -C 1. However, the value of the damping constant, 7,;, is proportional to k, and 
so the coupling strength controls how long it takes for the steady state to be reached. If the coupling between the 
SET and the resonator is switched off (i.e., n = 0), then Eqs. I|42|) and i|43|l break down because a steady-state for 
the resonator is never achieved as there is no damping. 

It is interesting to compare the expressions for the effective temperature and damping constant of the resonator 
with results for a quantum oscillator coupled to a single electrical tunnel junction (TJ)iiLSS Our results for a classical 
oscillator can be compared directly with calculations for a quantum oscillator when eVds ^ ^0 so that the electrons 
cause rapid dephasing of the resonator and induce the transition from quantum to classical behavior. 1 '' The effective 
temperature we obtain for a resonator coupled to a SET [Eq. (|44[)] is very similar to that obtained for a resonator 
coupled to a TJ, k^Tes — eVd s /2 (with the electron temperature set to zero in the electrodes). However, because 
(P)at(P)jv + i < 1/4, the heating induced by the SET will always be at least a factor of 1/2 less than that induced by 
the TJ. In fact, the effective temperatures obtained for both the TJ and SET can be commonly equated to one-half 
the ensemble-averaged energy dissipated by an electron due to tunnelling across a junction. For the TJ, the energy 
dissipated (for zero electron temperature in the electrodes) is always eVds- For the SET, on the other hand, the 
average energy dissipated by a tunnelling electron is Er(P)n + E l (P)n+i = 2E l (P)n+i — 2eVd s (-P) n (P) n+i for 

K < 1. 

The intrinsic damping constant for the TJ takes the form 7tj = /iC 2 /™, to lowest order in the coupling constant 
C, where /i is a constant which depends on the densities of states in the leads. In contrast, for the SET the intrinsic 
damping rate, 7; = (Re 3 /md 2 )(VgC g /vdsC-z) 2 , depends explicitly on the drain-source voltage. The difference in 
damping rates relates directly to the difference in electron tunnelling rates for a TJ and SET: for the TJ the electron 
tunnelling rates across the junction do not depend on the drain-source voltage, to lowest order, whereas for the SET 
the tunnelling rates always depend on the drain-source voltage in the orthodox regime. 

IV. RESONATOR WITH EXTRINSIC DAMPING AND AT FINITE TEMPERATURE 

Thus far, we have considered the nano-mechanical resonator and the SET as a closed system without including 
any external influences. Under these circumstances, and for sufficiently weak coupling to the resonator, the SET acts 
like a thermal reservoir which heats the oscillator and also gives rise to damping. The effective temperature in the 
steady-state is not a function of the coupling strength between the SET and the resonator, k, but the time taken to 
reach that steady-state does depend on the intrinsic damping rate which in turn depends on n. However, the resonator 
will have a finite quality- factor, even in the absence of the SET)24 and if the value of the associated damping rate, j e , 
is much larger than that due to the SET, then it will be essential to include these damping processes in the dynamics. 
Furthermore, both the resonator and the electrons will be at finite temperatures, set by their surroundings (which 
need not necessarily be the same). The effects of extrinsic damping and finite background temperature can be taken 
into account by modifying the master equation formalism we have developed so far. 
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A. Variances with Extrinsic Damping 

The equation of motion for the probability distribution is readily modified to include extrinsic damping of a harmonic 
oscillator arising from sources other than its interaction with the SET. We can obtain the correct form for the coupled 
master equations by adding the term 

Q 

— (j e uP N{N+1) ) 

to Eqs. I|16fl and (|17fl . where the extrinsic damping is modelled as a force on the oscillator of the form —j e u^ 

The mean-coordinate equations i|19l2UI21|) and the position of the fixed-point corresponding to the steady-state 
are unaffected by the extrinsic damping. However, this apparent lack of change belies a significant increase in the 
complexity of the underlying dynamics. With finite extrinsic damping, the steady-state values of (x) n and (x) n+i 
start to converge so that (x)n > and (x)n+i < xq(P) n+i- Furthermore, with extrinsic damping although the 
overall average velocity remains zero in the steady-state, the average velocities for given charge states are non-zero: 
(u) N = -(u) N+ i ^ 0. 

The variances of the steady-state distribution including extrinsic damping are obtained using the same method as 
before, and we find 

1 + & [(1-*)(1 + a) + e"] (45) 

= 1 + + (46) 

where a = j e Tt, and Sx 2 and Su 2 are the intrinsic variances (no extrinsic damping) given by Eqs. 142fl and 143fl . 
respectively. 

In the regime where k, e,a <C 1, the effect of the SET on the resonator will depend very sensitively on the ratio 
le/li — a./{t 2 n). When 7 e <C 7$ the extrinsic damping will not cause any serious modification of the resonator's 
behavior. In contrast, when 7 e ji the extrinsic damping will round-off the driving effect of the SET and we find, 

Sx 2 ~ Sx 2 — 



Su 2 ~ 5u 2 ^ 



and 



k B T eS ~ -eV ds (P) N (P) N+1 . (47) 

7e 



These expressions are consistent with the average back-action displacement noise derived in Ref. [2(j . 

The relatives value of the intrinsic and extrinsic damping constants are also important for understanding why a 
single- mode description of the resonator is sufficient. It turns out that the magnitude of the intrinsic damping constant 
is very much larger for the fundamental flexural mode of a resonator than for higher modes. It is for this reason that 
the behavior of the system is dominated by the fundamental mode and our single-mode treatment is appropriate, as 
we discuss in appendix C. 



B. Finite Oscillator Temperature 

The electrons in the SET and the resonator will be in contact with environmental degrees of freedom which will 
act as thermal reservoirs at finite temperature, but for the mesoscopic system we are considering, these temperatures 
will not necessarily be the same and their effect on the dynamics may be very different. The effect of the finite 
electron temperature, T e , will be to round out the tunnelling rates and thereby change the charge dynamics of the 
SET somewhat. However, from the full expression for electron tunnelling rate, Eq. QJ, we can see that this effect will 
be small so long as k^T e <C eVd s . 

The effect of a non-zero external bath temperature for the resonator, T r , can be gauged by adding in one further 
term to the coupled master equations, (|Td|) and lfT7|) . describing diffusion in velocity space. Thus in this case the extra 
terms take the form, 

9 , P ^ d ( G dP N(N+1) \ 
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so that without the coupling to the SET, the equation of motion for the resonator probability distribution is simply 
Kramer's equation^ For systems close to equilibrium, G is a constant related to the background temperature by 

G = 2mk B T rle . (48) 

In the regime where a, e, k <C 1 we can approximate the full equations for the variances by two the expressions 



6x 2 ~ — 



"0 



eV ds (P) N (P) N+1 , k B T r 



Su 2 ~ 1 
m 



!+7e/7i l+7i/7e 

eF ds (P) iV <P) i v +1 (l-/«) fc B i; 



l+7e/7i l+7i/7e 



(49) 
(50) 



The derivation of these relations rests on the use of the fluctuation-dissipation theorem [Eq. 148(1 ]. whose validity for 
the non-equilibrium SET-resonator system we have assumed without proof. Nevertheless, the resulting expressions 
for the variances have a very intuitive form. For k<1, the effective oscillator temperature satisfies 

T cS _ Tr_ + Tj_ /gjs 
Qef£ Qe Qi 

where the extrinsic, intrinsic, and effective quality factors are defined as Q e = uo/'Je, Qi = ^a/li — ^/( Ke )i and 
Q~g = Qe 1 + Q~ x , respectively. The intrinsic temperature, Ti, is given by Eq. Eq. fBTjl is just the steady 

state relation for a system in thermal contact with two reservoirs at temperatures T r and T^, and with the thermal 
conductivities of the pathways between the system and each reservoir proportional to the respective damping rates 
(and hence inversely proportional to the quality factors)^ 

The relative importance of heating due to the SET and that due to the resonator's environment is controlled by 
the ratios Qi/Q e and Ti/T e . When Qi -C Q e , we have T e g = Ti + T r Qi/Q e . As an example application, note that 
this equation predicts that it would be possible to actually cool the resonator to a temperature below T r by coupling 
it to an appropriate SET such that Tj < T r , 35 If, in addition to Qi <C Q e , we have that Tj ^> T r Qi/Q e , then the 
resonator hardly feels the effect of the external environment and its state is determined by the interaction with the 
SET electrons. In contrast, when Qi 3> Q e , we recover Eq. (|4Tjl provided T r <§; TiQ e /Qi. If instead T r ^> TiQ e /Qi, 
then the external environment dominates over the SET electrons and effectively determines the state of the resonator. 

In the regime where Qi 3> Qe-, but where the ratio Ti/T e is arbitrary, it is possible to describe the dynamics of the 
SET-resonator system using a very compact and intuitive formalism. Details of such an approach, which compliments 
that given here, will be described elsewhere^ 

V. CURRENT CHARACTERISTICS OF THE SET 

The analysis we have performed so far describes the effects of the electrons on the dynamics of the resonator. 
However, the electrons passing through the SET are of course influenced by the motion of the resonator and this is 
manifested as a change in the average current through the SET. 

In the absence of extrinsic damping, and at T r = T e = 0, the ensemble-averaged current through the SET can be 
calculated from the tunnelling rates and the steady-state probability distribution function: 

= el dx J du AE ^j x) P N+1 (x,u) (52) 

dx J du ^j^ P N {x,u) (53) 
= ~{P) n {P) n+x (1-k) (54) 

- ^ (55) 

Comparing this with the fixed gate result, Eq. (|HJ), we see that the resonator slightly reduces the current flowing 
through the SET. 

Including extrinsic damping, we find that the averaged current becomes 

a) = ah ( ? / 2(Q + K) T 1 Q(l tr ) ^ • (56) 



e 2 (a + k) + a(l + a)(l - k) 
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TABLE I: The resonant frequency, coupling strength and intrinsic quality factor for Si cantilevers with different lengths. 



1 in /xm 


1 


5 


10 


/,, in MHz 


240 


9.6 


2.4 


k/V* in V - * 


3 x 10~ Y 


9.8 x 10~ 4 


0.03 


QiV* in V" 


2.1 x 10' 


1.7 x 10 b 


2.1 x 10 4 



The average current is modified by external damping because it depends on the tunnelling rates through the individual 
junctions which themselves depend on the properties of the sub-ensembles through (x)jv and (x)jv+i- The external 
damping suppresses the motion of the resonator and this is reflected in the fact that for a 3> e 2 , kc 2 the average 
current is effectively reduced to that for the fixed-gate SET. 



VI. PRACTICAL CONSIDERATIONS 



The predicted effect of the classical back-action of the SET on the resonator using currently available technology 
depends crucially on the range of values of the parameters in our model which are accessible in practice. Typical, 
achievable values for the individual SET and resonator parameters are discussed in detail in appendix D. 

The maximum effective temperature of the resonator coupled to the SET, neglecting the effect of the extrinsic 
damping, is eVd s /&k&. Considering, for example, a typical total SET capacitance Cs = 0.5 fF and drain-source 
voltage Vds = 0.5e/Cs = 0.16 mV, the maximum effective temperature is about 0.5 K. Expressed in terms of the 
various parameters, the dimcnsionlcss coupling strength is 

\C S J V As muj z d 2 

Substituting in the expressions for C g , m, and luq appropriate for, say, Si and the above values for Cs and Vd s , Eq. I|57|) 
becomes 



w 



K = ixio- 12 v; 2 ^, (58) 

where the cantilever dimensions I, w, t, and cantilever-island gap d are expressed in micrometers, and the gate voltage 
V g is given in Volts. Similarly, the intrinsic quality factor is 

Q i = 1.3xl0 12 -^, (59) 

where again the dimensions are expressed in micrometers. Note the strong dependence on the cantilever length I and 
gap d. As might be expected, decreasing the gap and increasing the length strengthens the coupling between the 
SET and resonator and reduces the quality factor. This trend can be clearly seen for the range of cantilever length 
examples considered in Table [I] In each of the examples, the transverse cantilever dimensions are w = 0.25 /xm and 
t = 0.2 /im, and the gap is d = 0.1 /xm. 

For micron-sized resonators, the extrinsic quality factors, Q e = uJo/j ei of the fundamental modes typically lie in the 
ranged ~ 10 3 -10 4 , which, from the first column in Table [IJ will be much less than the quality factors associated with 
intrinsic damping even for gate voltages close to vacuum breakdown ~ 10 V (see appendix E). In this regime, for the 
one micron long example the effective temperature will be given, to a good approximation, by T e g = T r + TiQ e /Qi < 
T r + 0.24V 2 mK. Therefore, the maximum amount the SET can add to the cantilever temperature is about 24 mK for 
V g tuned to current peak maximum around the breakdown voltage = 10 V. The description of the SET back-action on 
the resonator in this regime coincides with the analysis given in Ref. 20. On the other hand, increasing the cantilever 
length by just a factor of ten, it is not difficult to be in the opposite, SET back-action dominated regime. With, e.g., 
V g = 5 V and Q e = 10 4 , we have Q e /Qi = 12 and T off = T 4 + 0.084T r ~ T, for T r < 12T;. 



VII. CONCLUSIONS AND DISCUSSION 



We have carried out a detailed analysis of the dynamics of a nano- mechanical resonator coupled to a SET. For the 
isolated SET-resonator system, the electrons tunnelling through the SET act like an effective thermal bath for the 
resonator. The variances of the stationary probability distribution for the resonator, and the dynamics associated 



14 



with its approach to the steady-state, allow us to assign an effective temperature and intrinsic damping constant due 
to the electrons, so long as the coupling between the resonator and the SET is small, i.e., k«1. 

The intrinsic damping constant controls the time taken for the isolated system to reach a steady-state. The effective 
temperature for the isolated SET-resonator system does not depend on the coupling between the resonator and the 
SET (for fixed values of (P)n and (P)jv+i), a fact which is somewhat surprising at first sight. However, although the 
average energy of the resonator in the steady-state is independent of the coupling, k, the time taken for the energy 
stored in the resonator to reach this average, constant level is determined by the coupling. In the limit that the 
SET-resonator coupling goes to zero, the time taken to reach a steady-state will diverge and so the whole concept of a 
steady-state for the coupled system and hence an effective temperature breaks down. The effects of extrinsic damping 
can be understood as trying to take energy out of the resonator, through an additional pathway, on a time-scale set 
by the extrinsic damping rate. Hence, the effective temperature of the resonator will now depend on the relative 
magnitudes of the extrinsic and intrinsic damping constants. 

We can also draw some conclusions from our results for the coupled dynamics which are more widely relevant to the 
study of nano-electromechanical systems. The idea that a mechanical resonator coupled to a mesoscopic conductor 
should have an intrinsic damping rate was demonstrated by Mozyrsky and Martinii and confirmed by Smirnov et 
alSiL The results we obtain for a classical resonator demonstrate the generality of this result. Thus far, many studies 
of mesoscopic conductors coupled to nano-mechanical resonators&ifiiiS do not discuss the possibility of damping of 
the resonator arising from its interaction with the electrons. Indeed, in many nano-electromechanical systems it may 
be that there is a non-zero intrinsic damping rate which affects the character of the steady-state. 

In this work, we have concentrated almost exclusively on how the properties of a nano-mechanical resonator are 
affected by coupling to a SET. Although we have calculated the effect of the coupling to the resonator on the average 
current through the SET, the spectrum of the current-noise has not yet been obtained, though work on this is in 
progress^ The current-noise spectrum will give important further information about the coupled dynamics of the 
system as it is very sensitive to correlations between the charge on the SET island and the position of the resonator. 

While preparing this work we became aware of Ref. |^} which considers the quantum dynamics of the SET-resonator 
system, but in a different regime to that considered here. 
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APPENDIX A: TOTAL ENERGY AND TUNNELLING RATES 

The circuit diagram for a nano-mechanical resonator coupled capacitively to the island of a SET is shown schemat- 
ically in Fig. ^ where a nano-mechanical resonator replaces the usual fixed voltage gate for the SET. In order to 
describe the dynamics of the coupled SET-resonator system, we require expressions for the total energy with N and 
N + 1 electrons on the SET island as well as for the tunnelling rates through the two junctions. 

1. Total energy 

The total energy of the coupled SET-resonator system for a given number of charges on the island is most easily 
obtained by generalizing the usual expression for the energy of charges on an SET island with a fixed voltage gate^i 
given in Sec. II. Here we treat the resonator as a single-mode harmonic oscillator. The generalization to several modes 
is discussed in appendix C. 

When the fixed gate is replaced by a harmonic oscillator, the gate capacitance, and hence the energy of the system, 
depends on the position of the oscillator. If the distance between the unperturbed oscillator and the SET island is 
d and the displacement, q, of the oscillator is always much less than this distance then the gate capacitance can be 
approximated by, 



C g (q)=C° g (l-q/d) 



(Al) 
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where C g is the capacitance of the un-displaced oscillator. Strictly-speaking we should write the approximation as 
Cg(q) = + qo(dC g /dqo) qo= o where go is the normal mode displacement of the fundamental mode, but we will 
neglect this detail as it leads to an unimportant geometrical correction. 

For nano-mechanical resonators, the values of C g are typically much less than the SET junction capacitances Cj. 
Hence, we will neglect the position dependence of Cs and deal only with position dependence in the term —2NN g . 
Thus, the charging energy of the SET island with N excess electrons can be written as, 

E ch = E c (V 2 - 2NN° g + 2iViV°!) , (A2) 

where E c — e 2 /(2Cs) and N® = —C^Vg/e (for simplicity we will drop the superscripts from now on as N g and C g 
will always take this form in what follows). 

The polarization charge induced by the gate can make the charging energies of states with N and N + 1 electrons 
on the SET island degenerate. For convenience, we define 

5N = N g -N-^. (A3) 

Using this notation, we find 

N 2 - 2NN g = SN 2 + i + 8N — N 2 (A4) 
(N + l) 2 - 2{N + l)N g = SN 2 + -- 5N - N 2 . (A5) 
Thus, we can write the total energy of the coupled SET-resonator system with N and N + 1 excess electrons as, 

H N = E C 6N + 2E c N g l ^N g - X - - 6Nj + J^- + ^mu 2 q 2 (A6) 

H N+1 = -E C 5N + 2£yV a | (n 3 + i - SN^j + ^ + \m^ 2 q 2 , (A7) 

respectively, where we have dropped terms which are constant in both. Making the change of coordinates x = q + Ax 
to eliminate the linear term in iJjv, we have 

\mulx 2 = \mujlq 2 + 2E c N g - J ( N g - SN - \ J + ^mu; 2 Ax 2 , (A8) 



2 



where 



Ax = 2EJ^ f _ 1 \ 

Expressing the total energies of the coupled system with N and N + 1 charges in terms of the new variable x, we 
obtain: 

v 2 1 

H N = E C SN + — + -muj 2 x 2 (A10) 
H N+1 = -E C 5N + \- -muj 2 (x - x ) 2 - -muj 2 xl (All) 

— III Zi 

-mujlxl(N g -5N - 1/2), 
where xq — —2E c N g /(mu> 2 .d) and we have again dropped constant terms. 



2. Tunnelling rates 

Within the orthodox model of single electron tunnelling, we can derive simple expressions for the tunnelling rates in 
either direction through the two SET junctions. Adopting the convention that electron tunnelling forwards corresponds 
to leftwards motion in Fig. ^ we can write down expressions for the four possible free energy differences: 

AE+ = -AEl = H N+1 -H N ~fi L (A12) 
AE+ = -AE R =H N - H N+1 + m , (A13) 
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where (J-l(r) 1S the chemical potential of the left(right)-hand lead. Using the expressions for Hm and ffjv+i obtained 
above we find, 

AE+ = -2E C SN - mw%x x - mu%xl(N g - 5N -1/2) + — ^ (A14) 
AE+ = 2E C SN + mulx x + mu>lxl{N g -8N -1/2) + -^, (A15) 
where we have used the fact that —^l — — eV^s/2. 

APPENDIX B: NUMERICAL INTEGRATION OF THE PROBABILITY DENSITY 

In this appendix we outline the way in which the coupled master equations, i|l(J|) and l|ll|) . can be integrated 
numerically. The coupled equations of motion for the probability distributions contain terms describing the free evo- 
lution of the resonator, which are deterministic, and terms describing the tunnelling of electrons which are essentially 
stochastic. The free evolution of the resonator is known exactly and the problem of obtaining the dynamics of the 
overall distribution is simplified considerably be disentangling this part of the motion from the electron tunnelling by 
transforming to an interaction picture. 

Formally, we can write the Liouville equation^! for Ppj(x,u;t), 

. dP N {x,u;t) 

i — = L N P N (x,u;t), (Bl) 

where Ln is the Liouville operator and a very similar equation for Pn+i{x, u; t). The dynamics of the resonator is 
simply that of a harmonic oscillator, albeit one whose equilibrium position shifts by Xq when an extra electron joins 
the island, together with a perturbation arising from tunnelling of electrons onto and off the island. Therefore, we 
can simplify the problem considerably by transforming into an interaction picture. We split the Liouville operator 
into two parts, 

L N = L%+ SL, (B2) 

where L t j^P^(x, u;t) = i {-Hat, Pn(x, u; t)} describes the free evolution of the oscillator with N charges on the island 
(for TV + 1 charges, we define L l j i+1 PN + i(x, u;t) — i {Hn + i, Pn + i(x, u; t)} by analogy) and SL describes the electron 

tunnelling processes. Defining Pjv(i) = e +lLNt PN(t), and a similar relation for N + 1, we have for the equations of 
motion in the interaction picture 

dPN{ 2 U]t) = ^{Q[AEi(x')}AEi(x') + e[AE^(x')]AE-(x')} (B3) 

xP N+1 (x + xq(1 — cos(wo*))j u + ujqXq sm(w t);t) 
1 

~Re 

xP/v(x, u; t) 

dP N+l {x,u-i) = _^ {e[AS + ( ^ )]Ai? + ( ^ ) + e[Ai? - ( ^ )]A ^ ( ^ )} (B5) 
xP N+1 (x,u;t) 

+ {e[AE+(x")}AEt(x") + e[AE R (x")]AE R (x")} (B6) 



2 {Q[AE+(x')]AE+(x') + Q[AE- L {x')]AEl{x')} (B4) 



where 



Re 2 

xPn(x — xo(l — cos(woi)), u — luqXo sin(o;ot); t) 



u 

x cos(a->oi) sin(o;ot) 

u 

xq + (x — Xq) cos(u>ot) sin(woi) 

UJf) 



(B7) 

(B8) 
(B9) 



17 



Thus the transformation to the interaction picture eliminates the free-evolution of the oscillator and leaves us with a 
pair of simple rate equations which can be integrated numerically in a straightforward way. 

Integrating these equations numerically on a finite grid and then performing the reverse transformation back from 
the interaction picture to obtain the full time dependence, we can calculate the probability distribution as a function 
of time. The accuracy of the calculation is set primarily by the size of the grid used, but convergence of the results 
can be confirmed by performing the calculation on grids with progressively smaller mesh sizes. 



In the main body of this paper we treat the resonator as a single-mode harmonic oscillator and neglect the effect 
of all but the fundamental mode. In this appendix we show why such an approximate description works well. 

In most studies of nano-electromechanical systems in which the mechanical component consists of a beam or 
cantilever, which has a spectrum of vibrational modes, all but the fundamental flexural mode is usually neglected. For 
systems where the coupling between modes is weak, this is a sensible approximation as it turns out that geometrical 
factors mean that the coupling between a normal mode of a beam or cantilever decreases rapidly with increasing 
frequency^ 

In the system we consider here, the situation is a little different as the electrons passing through the SET act as a 
thermal bath and heat the resonator up to an effective temperature that is independent of the coupling. Therefore, in 
the absence of extrinsic damping each of the normal modes would have the same effective temperature. The overall 
variance in the position of the resonator would then be obtained by summing over the contribution from each of the 
normal modes. However, when we include the effect of external damping (discussed in Sec. IV,) it is found that the 
effective temperature for a single mode depends on the ratio of intrinsic to extrinsic damping rates: 7i/7 e - When the 
ratio 7i/7e is <C 1, the extrinsic damping dominates and the SET has little effect on the mode. However, for ji/j e < 1 
the extrinsic damping does not fully counteract the thermal-like excitation of the mode due to the SET. 

The intrinsic damping, 7^, can be written in terms of the basic parameters of the system as 



This formula follows from the approximation we made in appendix A where we neglected the distinction between the 
position of the center of mass of the resonator, q, and the normal-mode displacement of the fundamental mode. In 
fact, the linear correction to the gate capacitance should be written as {dC g /dqo) qo= Q = —xCg/d, where go is the 
fundamental mode position co-ordinate and x is a geometrical correction factor which we have neglected as it is of 
order unity for the fundamental mode. 

In general, there will be an intrinsic damping constant for each normal mode of the resonator, n, of the form 



where q n is the displacement of the n-th normal mode, 15,42 However, the intrinsic damping constant for higher modes 
will always be much less than that of the fundamental mode because of their shorter wavelengths and the presence of 
nodes, which significantly reduce the value of the coupling, given by (dC g /dq n ) qn= Q.' 12 Furthermore, there is reason 
to expect that the extrinsic damping, 7™, will be greater for higher- frequency modes^ leading to a further reduction 
in the ratio 7™/7™ , and hence the effect of these modes on the overall motion of the resonator. 

The current characteristics of the SET and the associated noise spectrum are also much more strongly affected by 
the fundamental mode of the resonator than by higher modes. This is again because of the weak-coupling between 
the resonator and the SET. 

It should be noted that the treatment here could be extended to include the effect of higher modes of the resonator. 
The calculations presented in the main text consider a single mode, and the correct multi-mode results could easily be 
obtained by adding up the contributions of each of the modes, although the appropriate coupling constant, k (which 
in general would depend on the frequency and geometry of the mode), would need to be used for each of the modes. 
The only assumption in our approach is that the normal modes of the resonator are effectively uncoupled from each 
other on all the time-scales of interest. 



APPENDIX C: HIGHER MODES OF THE RESONATOR 




(CI) 




(C2) 



APPENDIX D: PARAMETER VALUES 



In this appendix we discuss the physical values of the various parameters in our model of the SET-resonator system 
which could be achieved using available technology. We concentrate on the parameter values which would be necessary 
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for the system to operate in the regime we have explored theoretically, where the SET is well-described by the orthodox 
model, and the drain-source voltage is large enough that the resonator behaves classically. 

For the orthodox model to apply, the charging energy, E c = e 2 /(2Cs), must dominate over the electron thermal 
energy scale, k^T e . The total capacitance of the SET, Cs, is dominated by the capacitance of the tunnel junctions, 
Cj, typically much larger than that of the mechanical gate, C g . A typical choice of Cj — 0.25 fF gives E c /k-Q — 2 K so 
that at electron temperatures of a few tens of mK, the charging energy clearly dominates. The junction resistances, 
R, of the SET must be set relatively high, e.g., R — 100 kS7, in order to suppress higher-order electron tunnelling 
processes which cannot be described within the orthodox model. The drain-source voltage is typically chosen to be 
of order the charging energy eVd s ~ e 2 /(2Cs), as it is in this regime that the SET conductance is most sensitive to 
changes in the gate capacitance. Thus for Cj = 0.25fF, we have, e.g., Vd s = e/(2Cs) = 0.16 mV and the electron 
tunnelling time, r t , is consequently 2RC% = 10~ 10 s. 

The properties of the nano-mechanical resonator used as a voltage-gate can be varied to quite a high degree. The 
frequency of a nano-mechanical resonator is determined by the material it is made from and its geometrical form. 
For a cantilever, a beam which is clamped at one end and free at the other, the fundamental flexural mode frequency 
is given by ojq = 1.02(i// 2 )^/ Ej p, where I is the length, t the thickness, E is Young's modulus, and p the mass 
density of the material from which the resonator is fabricated. Resonators fabricated from semiconductors typically 
have y/ E/p in the ranged ~ 10 3 -10 4 ms" 1 . For the example of a Si cantilever with E = 1.33 x 10 11 Nm" 2 and 
p = 2.33 x 10 3 kgm -3 , the fundamental flexural mode frequency is /o = 1.2t/l 2 GHz, with the dimensions expressed 
in micrometers. The gate capacitance of a nano-mechanical resonator is also determined by geometry. For a vacuum 
gap, C g = Xeolw/d, where w is the cantilever width, d the cantilever-SET island gap, and A is a geometrical factor 
(which accounts, e.g., for the fact that the metal electrode may have smaller area than the cantilever surface Iw). 
Considering, e.g., a resonator with dimensions I = 1 /im, w = 0.25 psm, t = 0.2 /im, d = 0.1 /im, and A = 1/3, we 
have C g — 74 aF and /o = 244 MHz. This corresponds to e = 0.02 for r t = 10~ 10 s. The resonator would be well 
within the classical regime as the energy quanta for this frequency are hujQ = 0.15/ieV, very much less than the energy 
associated with the above drain-source voltage. 

The magnitude of the gate voltage, V g , is limited by the need to avoid electrostatic breakdown in the dielectric 
between the SET and the resonator ; 14 ' 19 ! 44 for vacuum the breakdown voltage is about 10 8 V/m^ translating into a 
maximum voltage of about 10 V for a resonator-SET island separation d = 0.1 /im. 
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